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Ph 1 Abstract 

-P .... 

We consider an application of ideas of inertial confinement from plasma physics for the gas of 

charged particles interacting at high energy. We formulate this problem as a problem of equilib- 
rium state for the "internal" (slow) particles in some additional potential field created by "external" 
f-^. ■ (fast) particles. In this formulation we find an equation of state for the "internal" particles. We 

also consider a shock wave description of the process discussing a self-similar solution of the Eu- 
ler's equations for the fields of interest and we solve Boltzmann equation for self-consistent field 



(Vlasov's equations) in linear approximation determining an electrical properties of collisionless 



plasma. Finally we discuss the possibility of the applications of the methods of the approach for 
the investigation of the early stages evolution of quark-gluon plasma. 
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1 Introduction 

A collision of heavy nuclei at LHC offers new possibilities to explore a new state of matter, the Quark 
Gluon Plasma state, which is a dense state of strongly interacting quarks and gluons. Initial stage of 
scattering of two nuclei is an initial condition for further evolution of the bulk of quarks and gluons in 
the form of the ideal fluid, see [1] . From this point of view it is interesting to understand a creation 
of the very dense areas of the mater at initial stages of scattering. 

The process of nuclei-nuclei scattering at high energy is a highly non-equilibrium process due to 
the presence of the partons with different distribution over the particles velocities (rapidities), see [2] 
for example. In order to describe processes of interactions between these partons the approaches of 
[3] and [1] can be used. We assume that part of the bulk of the particles form drops of the partons 
of very high density due to the first initial nuclei fronts interactions. Approximately half of the total 
initial energy of the interacting protons or nuclei is going there. These particles we can call as internal 
("slow") particles: they loose their velocities during the initial processes of interactions at collision. 
This picture is due to [1] , at initial stage of interactions the " slow" particles are in the rest in central 
mass frame, further expansion of these drops is the matter of the hydrodynamical description. There 
are also partons which we can call as an external or "fast" ones, they remain relativistic, i.e. they 
almost do not loose kinetic energy and they do not change their quantum numbers. The influence 
of these particles in the process of the hot drops compression we can account by their action on the 
"slow" partons. 

The separation between the internal and external particles, therefore, is based on their sequence in 
the involvement in the process of interactions and their relative velocities in the laboratory frame. The 
kinetic energy of the "slow" particles is heating the volume of the drop whereas the "fast" particles 
almost preserve their velocities. Thereby the main difference between these "slow" and "fast" particles 
comes from the value of their relative rapidity after the collision - the secondary " fast" particles have 
a large rapidity in relation to the bulk of stopped "slow" particles immediately after the collision of 
target and projectile. Now, following to [3], the interactions of the internal particles with external ones 
we can describe as an equilibrium process for the internal particles in some additional potential field. 
Namely, the influence of the "fast" particles we can reduce to some effective potential field included 
in the description of the stopped matter state. The advantage of this approach is clear, in this way 
we reduce the non-equilibrium process to the equilibrium one with the possibility to apply a wide 
variety of approaches knowing for equilibrium description of the matter states. The whole process 
description, therefore, is similar to the effect of inertial confinement in plasma physics, see [5], where 
for the creation of the requested density of the matter needed for the thermonuclear fusion a radiation 
pressure on an external area of a target is used, see the Fig(TJ 

In our paper, in order to understand a qualitative behavior of the internal particles, we consider 
an oversimplified model of interacting particle. We assume that the gas of our particles is the gas of 
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Figure 1: Inertial confinement in the thermonuclear plasma fusion. 



charged disks with some additional inter-particle interactions which depend on rapidity variable. This 
rapidity variable characterize the secondary (" fast" ) particles and related with the total energy of the 
process. The choose of this variable as the main variable which separate stopped partons from the 
fast secondaries in the process is dictated by the QCD Pomeron approach in the high energy physics, 
see [6]. Indeed, the process of interactions between internal and external particles is very similar to 
the high energy interactions governed by the Pomeron amplitude. In some sense the processes are the 
same excluding the Pomeron sources terms of the processes. For usual high energy interactions these 
sources are the target and projectile, whereas in our approximation the sources are the flow of the 
fast partons and drop of the stopped partons. Nevertheless, in both cases the main variable of the 
external potential must be a rapidity variable. The value of the potential, therefore, is proportional 
to the total energy of the process what is also a physically clear and acceptable assumption. 

The main aims of our consideration, which can be useful for the better understanding of the 
properties of quark-gluon plasma, are the understanding of qualitative behavior of matter at initial 
stage of interactions. In particular we will discuss the transition from gas to liquid description of the 
interacting matter and some transport properties of the system which can be relevant also at the case of 
real quark-gluon plasma. The Section 1 is an introduction of the paper. In Section 2 of the manuscript 
we consider the interactions between our "toy" particles with very general pairwise potential in two 
dimensional space. The additional external field is chosen to be depended on the rapidity variable 
and it is the only characteristic of the external particles in the model. For our purposes the form of 
the potential is unimportant and we chose some simple one. The following subsections of this section 
are dedicated to the Van der Vaals equation of state of the gas. Section 3 describes thermodynamical 
properties of the gas. In particular, we discuss transport characteristic of the gas basing on Enskog 



approach, see [7]. Section 4 is dedicated to the shock wave description of the process of interactions 
of external particles with internal ones. Namely, the process of external particles influence on the 
internal ones is described as a convergence to the center of area of the internal particles shock wave. 
We consider Euler's equations for the case of this shock wave and discuss the obtained results of this 
approach following mainly to the results of paper [8]. In Section 5 we generalize the problem and 



consider Boltzmann equation in Vlasov's 
describe an equilibrium state of the system 



brmulation for collisonless plasma in three dimensions and 
Last section is dedicated to the discussion of the obtained 



i 



results, there the conclusion of the paper is given. 

2 Equation of the state of the gas 

2.1 The potential of pairwise interactions 

In order to write the equation of state of the gas we need to determine the pairwise potential between 
two particles, see Appendix A. So, a potential energy of our system we define as a sum of the two 
parts: 

U(b) = J f(b-x)V l (b)f{x)d 2 x + V 2 {b). (1) 

Here V\(b) is the usual pairwise potential arising between the particles at the mutual distance b: 

V 1 (b 12 ) = -A/bu + B/b? 2 (2) 

with a > 1 (in our further calculations we take a = 2) and A, B are some constants.. The function 
f(x) is a distribution function of the particles which we take equal 



2 /p2 



P -x 2 /H 

m = -^r. (3) 

see |10j for justification of that choose. The potential energy term V 2 (b) , absent in Eq. (|A.ip . describes 
the potential energy of a single particle inside of some region (with radius R) in the field created by the 
particles outside the region, i.e. this term accounts an influence of the fast, energetic particles outside 
of the region with radius R on the particles inside the regioio Namely, following the proposition of [3] , 
we reduce a non-equilibrium process of interactions to the equilibrium interactions of only part of the 
particles in the potential field of external force created by another part of the particles. Following our 
assumption about the Pomeron like structure of this external field, an additional potential field V2(b) 
depends, therefore, only on the position (impact parameter) b of the particle inside of the region and 
rapidity Y = In ( j- J of outside particles in respect to stopped, "slow" , ones, where s as a squared total 



1 In Appendix B, related with the Section 5, we investigate a deviation from the equilibrium state as well 

2 The value of R is somehow arbitrary, it can be considered as radius of proton for example, or as a characteristic 

value of some hotspot inside the collision region. 



energy of the "fast" particles in c.m.f.. This field creates additional pressure on the "inside" particles 
and could be considered as an analogue of the inertial confinement effect in physics of plasma. For 
simplicity we choose the form for this term similar to what we have in the phenomenological Pomeron 
approaches, see [TU] : 

V2(b) = -j^e Y 8(R 2 -b 2 ) (4) 

where 9 is the step function and C is a constant. In general, the separation to the "internal" and 
"external" particles depends on a time evolution of the stages of collision, we assume that some 
region of "slow" particles with high density is produced due very initial stage of nuclei interactions. 
These particles inside the region we call as " internal" , whereas " fast" secondary particles we can call 
as "external" ones. If, due the interactions, these "fast" particles transfer all the momenta to the 
particles inside the drop and they also stop then anyway we consider them as external to for the 
particles inside the initial drop. 

For the potential energy we obtain finally the following expression: 

p -b 2 /2R 2 p -b 2 /2R 2 n 

U(b) = e -^- w V l {b) + V 2 (b) = e -^- w (-A/b + B/b 2 ) - ^ e Y 6{R 2 -b 2 ). (5) 

Considering an integral Eq. (|A.12p with the potential Eq. ©, we see that we can approximate the 
integral by the following expression 

/rbo poo 

d 2 b{exp(-pU l2 ) - 1) « - / d 2 b-P U 12 d 2 b, (6) 

JO Jb 

where bo is the value of b where 

U 12 (b o )^0. (7) 

Thereby, the value of &o i n Eq. (|A.12j) we find as roots of the following simple equation: 

" J + | " ° eY = ° ^ 

which gives in the limit of large Y (high energy limit): 

bo(Y) « e~ Y ' 2 U%~± e- Y l 2 \ = b (Y) e^ 2 . (9) 

Clearly, at very large values of final rapidity Y we obtain that bo <C B/A, that means bo <§i 1. 
The reason of this decrease of the inter-particle distance is the presence of the additional pressure 
introduced in the model. 

2.2 Equation of state of the gas 

Now, we ready to calculate the integral Eq. (|A.12p : 

J 2 = f d 2 b (exp {-p U 12 ) - 1 ) w - / ° d 2 b - p [ U{b) d 2 b = (10) 

J JO Jb 

5 



Performing the integration and keeping in the integrals only leading terms of answer's expansion at 
large Y (bo(Y) € 1) we obtain: 



f3e Y C 

+ ^2~r 2 ~ 

Using shorter notations we write the answer as 



(R 2 -bl(Y)). 



J 2 « -nb 2 + $ (l(Y) + B(Y) + C{Y)e Y ) . (12) 

Therefore, the state equation of the gas in our case has the following form: 

P =^f-^- {-« b o(Y) + P {MY) + B(Y) + C{Y)e Y )) , (13) 

or, using the shorter notations again, we have: 

P = ^f~^f (-7T&8QO + Pa(Y)) . (14) 

Introducing the gas density variable 

N . , 

P = Jl m -n. ( 15 ) 

we rewrite the state equation in the following form: 

(P + AOO) (l - pnb 2 (Y)) = pT . (16) 

This equation of state of the gas is a modified Van der Vaals equation with coefficients of equation 
which are depend on the rapidity of the process. 

2.3 Critical parameters of the equation of state and thermodynamical parameters 
of the system 

The critical parameters of the equation of the state we find similarly to the usual definition of critical 
parameters in the Van der Vaals equation 

3 , P "ct'KUq +T cr P cr 3 

~ P + —& ~ P U + T2 = (Per ~ P) , ( 17 ) 

7r Og an 6q an Oq 

that gives the well-known answers 

Pcr ~ 3nb 2 (Y) ' nr ~ 27 tt bl(Y) ' ±cr ~ 27 n b 2 (Y) ' [Lb) 

Now, using Eq. ([9]) and approximating a(Y) ~ Ce we obtain 

o ~ -^—e Y P - ° e 5Y/2 T - 8 ° e 5Y/2 (19) 

3nbl 21nb\ 21 n b\ 

We see that all the values of the critical parameters in the critical point of gas-liquid phase transition 

are much higher then in the case of usual Van der Vaals equation. 



3 Entropy and transport properties of the gas 

The state equation of the gas, Eq. (|16p . we rewrite it in the following form 

P = Z(Y,T)pT, (20) 

where 

Z(Y,T) = l —^— - /3pa(Y) (21) 

1 - pirb£(Y) 

The internal energy (per particle in one mole) of the gas with two transverse degrees of freedom , 

therefore, is given by usual caloric equation 

U = PS = P (22) 

N Z(Y,T) pZ(Y,T) v ; 

where the Gibbs equation for U takes place 

P 

dU = Tds +^dp (23) 

P 2 



with s as an entropy. Using the caloric equation Eq. (|22p together with the state equation Eq. (|20p . 
we could rewrite the Gibbs equation for the entropy Eq. (|23p 



ds = -dU ^dp (24) 

T Tp 2 H y ' 

in the following form: 

dP dr/ 1 dZ dri 

ds = — - — dr) - Z— (25) 

P rj Z drj r\ 

where m is a particle mass and rj(Y) = irb^Y) p we could call as a packing factor of the problem. 

Integrating on of the Gibbs relation Eq. (f25j) gives for the entropy per particle in the gas 

P 

s = In j — YT\ \" + s o ! (26) 

7]Z(r]) expi^f -^-dnj 

see |S]. At high energy limit an integrating over r] gives for the entropy 

s = ln P(1 7 T?) + s . (27) 

rj z 

At high energy limit, when t/<l, using Eq. ([9]) we obtain the following expression for the entropy 

s = In— + 2Y - 2ln(7rb 2 {Y)p) + s , (28) 

which linearly depends on the total rapidity Y of the process. 

The equation of state, Eq. (|16|) or Eq. (|2Up . we can rewrite in the following form at the limit of 
high energy: 

P = pT (l + rjZ(rj)) , (29) 



where 

Z(v) = -r^— ■ (30) 

1 — T) 

In this equation the compressibility factor Z{rj) from Eq. (j29|) coincides with the Enskog factor in 
Eq. (|20p . Therefore, the transport properties of the gas, which are characterizing by shear viscosity 
C, bulk viscosity £ and heat conductivity k we take in Enskog approximation from [7J: 

c = Cor/ (^y + 1 + - 87??z(r?) ) ' (31) 

£ = £ r? (1.257^(77)) , (32) 

K = K0?? (^) + 5 + °- 87r l Z (v)) ■ (33) 

The Enskog approximation, therefore, allows to find the relative changes of the transport coefficients 
of the system under the pressure whereas the initial values of these coefficients are given. At the limit 
of high energy, whereas rj <C 1, keeping only leading in rj terms, we obtain: 

C « Co , (34) 

£ = 1.25 to V 2 Z(r)) « 1.25£ e" 2y (vr^(y) p ) 2 (i + e - y ff ^(y) p) , (35) 

k = K (l + ~r)J . (36) 

We see, that the bulk viscosity £ is small at the limit of large final rapidity Y in comparison to the 
case when Y = for our gas. 

Now we can consider the obtained results for the entropy and transport coefficients properties from 
the point of view of results for quark-gluon plasma models. In [9] was obtained that £ /(, ratio has 
a maximum at some critical temperatures T c . Considering the same ratio constructed from Eq. ()34|) - 
Eq. (|35p values we obtain that this ratio is equal to 

| = A L (37) 

and that it is small. Moreover, considering the critical parameters of gas-liquid transition Eq. (|18p we 
see, that this ratio has a maximum at p = p cr and at critical temperature Tjf|: 

* - L25e ° f38l 

c"^cT' (38) 

which is similarly to the properties of the same ratio in the quark-gluon plasma model of [9]. Thereby 
our simple model demonstrate properties which are similar to the properties of much more complicated 
models of quark-gluon plasma. Physically it means, that in Enskog approximation the bulk viscosity 



3 The critical temperatures in [9] and here have different meanings, but both of them define higher temperature 
comparing to initial one. 



is proportional to the density of the state and the maximum values of density and bulk viscosity are 
achieved at liquid state of the system. Further, when the state of non-compressible liquid is achieved 
due the evolution, the density and bulk viscosity of the liquid will remain constants. 

Additionally it is interesting to note that whereas in our model the entropy is a growing function 
of the total rapidity of the process, see Eq. ([25]) . the shear viscosity Eq. ([3~4"]) is a constant. It provides 
a small Q / s ratio 

C Co , QfV . 

- oc - (39) 

as it found in hydrodynamical description of quark-gluon plasma evolution, see in [T] for example. 



4 Shock wave point of view 

In this section we proceed following to main results of paper [8J. Our flow of compressed interacting 
particles could be described by the following Euler's equations for the radial velocity, mass density 

d t a + dJau) + - pu =0, (40) 

r 

dtu + ud r u H d r a = , (41) 

a 

and, assuming an absence of dissipative processes for our gas, for the entropy 

d t s + ud r s = 0. (42) 

Here u is the radial velocity of the flow, a = p m is the mass density of the gas and s is the entropy 
from Eq. (|27jl : 

with 

C(a) = - ln ( 1 - a n b 2 (Y)/rnj 2 . (44) 

There are three partial differential equations for the main fields which characterize the flow of in- 
teracting particles under the external pressure. Solution of this system is considered in the next 
subsection. 

4.1 Shock wave description of the process: self-similar solution 

In this section we consider the process of gas compression as the process of propagating convergent 
radial shock wave, which can be described in term of scaling variable 

r shock A (t - t) a 



where A and a are some constants to be determined further, an implosion of the shock wave occurs 
at time t = to which corresponds to £ = oo , and the front of shock wave is described by the radius: 

r shock = A (t - t) a . (46) 

The fields of interests are assumed to have the following form 

a = a G(O, (47) 

« = t^t no , (48) 

t — to 

with functions G(£) V (£) W(£) to be found from the equations Eq. (00]) - Eq. (pl3|) which we rewrite in 
the following form: 

(,V + (V-l) d -^ + 2V = 0, (50) 



dint: v ' dln£ 

V ' dlntl 2 dlntl 2 dint: V« / 

dW dlnG dC dG 2 1-aV , , 

+ --T^^ = 0- (52) 



dln£ dln£ dGdlnf a 1-F 
The system of differential equations Eq. (|50p - Eq. ()52[) can be linearized by expansion of the functions 
V, W around V(£ — > oo) — > and W(£, — > oo) — > that corresponds to t —> to time in Eq. (]4"8]) - 
Eq. (|49|) , see more details in [8|. The solutions of the linearized equations can be easily found, they 
are 

"«^. («) 

with constants ify and .KV must be found from the full equations. Concerning the G function from 
Eq. (|47p we note, that this function remains finite at £ — > oo: 

dlnG 



(55) 



dint: 

when V(£ ->• oo) ->• and M^(£ — ► oo) — > 0. 

In order to solve the equations Eq. ([50]) - Eq . ([52]) . we need initial values of our functions V(£ = 
1), Z(£ = 1), G(£ = 1) which could be determined from the matching equations on the discontinuities 
of the flow: 

<7i u\ = a 2 u 2 , (56) 

Pi + a 1 u\ = P 2 + o- 2 u§ , (57) 

P] mv% Po rauk , . 

Pi 2 p 2 2 

10 



with subscript 1 for the quantities before the shock and subscript 2 for the quantities after. The 
solutions of these equations in the limit Pi S> -Pi were found in |8] and they are the following: 

" - 2aiU * (60) 



2 + Z( m ) 
1u\ 

U2 - Ul = -YT^m)> (61) 

Q"" shock /rri\ 

U\ = -T Shock = — • (62) 

t — to 

Additional equation which relates the densities of the flow is Eq. ([17)1 : 

a 2 = a G(l) = aiG(l) (63) 

that together with Eq. ([5T)j) gives: 

G(1) = » + Z {m J o( Y)G(l)) = ' + 2 (' - »"8«°W) ' < 64 > 

Therefore, we obtain for G(l): 

nn\ = 

1 + p 7Tb 2 (Y) 

that in the high energy limit F» 1 gives poirb^Y) <C 1 and 



G(l) = 7— TI^T" (65) 



G(l) « 3 (66) 

for any value of po . The functions V(l) and W(l) could be found as well and they have the following 
form: 

V(l) = 1 - jjLj (67) 



and 



^(1) = 2^ I ^- (68) 



A numerical values of a = 0.8 and G(oo) = 4.6 are calculated in [8]. 

4.2 Parameters of the solution and fluid state of the gas 

In order to calculate parameter A from the Eq. (|46p we use the condition on the front of shock wave 
at initial time t = : 

r(t = 0) = At$ = R (69) 

with R from Eq. ([3]). So far we obtain 

•A = £- (70) 



'0 
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On another hand, the velocity of the shock wave at the initial moment of the time may be found from 

Eq. gED: 

c sh = V(l) ^ , (71) 

to 

that gives 

t = aV(l) — . (72) 

Csh 



R ( c sh " \ r) ] - , , 



Finally we obtain: 

A = is = {^m )«'-"■ (73 » 

Now we consider the pressure, Eq. (|49p . achieved in the system of interest with known solution for 
G,V,W functions: 

P - W^W° G ® W ® = 0? " < f iKw A'"r'-'" (74) 

Before the compression of the gas the state equation of the gas is the usual one: 

To = ^ , (75) 

whereas at time of compression Eq. (|20p holds: 

TiQ= l&- (76) 

Therefore, with the pressure from Eq. (|74j) we obtain for the temperature ratio: 

f^TOM GK) ^ 2/ ^ 2/ "- (77) 



Simplifying the expression we obtain finally 



mc~ 



R\ 2 '"- 2 



The maximum temperature in the center of area is achieved when £ — >■ oo at r — Y 0. Taking this 
limit and changing r in the expression by the minimal possible distance Eq. ([9]), that smooth out the 
divergence, we obtain the maximum temperature achieved in the system in comparison to the initial 
temperature Tq = Pq/pq: 



r 2 f p \2/«-2 

T(e = oo) =-. e^ — "' C '* „ T „^, K w (tJ±A (79) 



2Z(p 7r62(y )G( oo))y(l)2 ^00/ 



with 

1 — a . 

A = « 0.25. (80) 

a 

From another hand, we could compare the maximum achieved temperature with the temperature 
To = T(£ = 1) on the edge of shock wave. In this case we have: 

T(fl _ Zfa(l)) ^V /tt ~ 2 (81) 



T(l) G(l)Z(r?(0) Vr 
12 



and again smoothing out the divergence we obtain: 

r(oo) _ ay Z(p irbt(Y)G(l)) ( R \ 2/ °^ 
T(l) ' G(l)Z(p Q nb*(Y)G(oo)) \b (Y) ) ' { ' 

The resulting expressions Eq. (|79p . Eq. (|82p is interesting from the following point of view. Taking 
ratio of temperatures in Eq. (|82p from Eq. (|19|) as 

(83) 





T(oo) 
T(l) 


oc 


T 

± cr 




KS 


p 5Y/2 




T(Y = 


0) 




arrive 


to the following expression: 
















R 


-~r* 


Y 


> 


1 



(84) 
bo(Y) 

with 



It is well known, see [TT] for example, that when the characteristic parameter of the system, which 
is R in our case, is larger than the average distance between particle, which is bo(Y) in our case 
approximately, then we can consider a fluid flow instead initial gas statqj. Thereby we obtain an 
important analytical result: at large enough external pressure, which is characterized by the rapidity 
in the present model, we can trace a dependence of the transition from gas phase to the liquid phase 
on the given external parameters. This is interesting especially because, as it wide accepted, the initial 
stages of quark-gluon plasma evolution is well described in the framework of hydrodynamical approach 
, for the review of the subject see [1] and [5j for the description of hydrodynamical expansion of the 
state. 

5 Boltzmann equation for self-consistent field 

In this section we consider a kinetic approach to our problem, which we reformulate in comparison to 
the Section 2 - Section 4. Namely, now we consider a gas of charged particles of the same kind interact- 
ing now in 3 dimensional spaced which are under the pressure of external particles. When dissipative 
processes are negligible and local equilibrium is achieved, the Vlasov approach to the Boltzmann 
equation is valid. We assume that this is the case and, thereby, we consider the following Boltzmann 
equation for the one-particle distribution function in Vlasov's approximation of self-consistent fields: 

dJh^A +v d J^A + ^^1 = (86) 

ot or op 



4 The ratio opposite to Eq. (|84[1 ratio is an analog of Knudsen number in hydrodynamics. 

5 We formulate the problem in 3 dimensions because of importance of longitudinal dimension in real high energy 

interactions. The radial symmetry of the problem in this formulation is preserved as well and in the following we assume 

that all vectors of interest are radial. Therefore, the vector notations will be use only when it will be need in. 

13 



where we adopted the radial symmetry of the problem and where 

F(r) __£«w + , B __a«W »£W 

or or or 

is a force which consists of two terms arose from self-consistent electrical potential and external 
potential is similar to the one from Eq. Q: 

V 2 (r) = JjLpe Y f(r). (88) 

Here /(r) is some function of r and E = — grad(f>o(r) in Eq. (|87p is a self-consistent electric field 
created by our charged particles. This field is assumed to be weak enough in order to justify a linear 
approximation for the distribution function. We assume also that the magnetic field is small and we 
neglect the magnetic filed in our calculations. The Vlasov's system of equations include Maxwell's 
equations for the electric field: 

rotE = 0, div E = A-Kqu / f(r,p,t)d p, 

where n is a particle's density. 

5.1 Linear approximation: equilibrium state 

We solve the Boltzmann equation Eq. (|86j) by approximation of distribution function with some linear 
correction to the equilibrium state 

f(r,p,t) = fo(r,p) + fi(r,p,t), h « f (90) 

and some linear correction to the constant electrical field 

E = E° + E 1 . (91) 

We are looking for an equilibriums state at the presence of external pressure in initial condition for the 

equation. Therefore, the Boltzmann equation which describes this equilibrium state has the following 

form: 

^ + F(r) !M = 0> (92) 

drop 

whose solution is the Boltzmann-Maxwell distribution function: 

fo(r,p) = Tj^e 2mk BT k B T k B T _ ( 93 \ 

(2irk B Ty 12 



6 In general, this state is quasi-equilibrium, the external pressure acts during finite period of time. Nevertheless, we 
assume, that this time period is long enough and we could consider this state as an equilibrium one during the early 
stages of system's evolution which we are interesting in. 
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The Maxwell equation for the electric field Eq. (|89p . therefore, reduces to the Poisson equation which 
has the following form in the limit of high temperatures (large kinetic energies of the particles in 
comparison to the potential energy): 



or 



■i . >M r ) _ A ( -, v 2(r) 



A(f)o(r) - 4irq 2 n-^ = -inqn f 1 - -^ ) . (95) 



We could rewrite it in the following form: 



A^ (r)-^ = -Qo + ^, (96) 

r D Q r D 



with Qq = w J n as the charge density and 



2 k B T 



r D = , o (97) 

as the Debye length. Taking, as simplest example, the potential Eq. ([1]) in form of Eq. (JU) 

V 2 (r) = ^e Y 9(R 2 -r 2 ), (98) 

we are looking for two solutions in two different regions: 

1. Region where r > R. In this region the potential Eq. (pQ) is zero and we obtain the usual 
screening equation for the electric potential: 

A0 o (r) - ^P- + Q = 0, (99) 



whose solution is 



-r I r D 

(r) = Qor 2 D + C . (100) 



2. Region where r < R. In this region we have 



AMr) -^ + Qo~ ~ ° 2 e Y = (101) 



' D *"" ^'D 

with solution for the potential 



— r I rj3 n 

Mr) = Qor 2 D + C' — ^^7 e Y . (102) 



Matching the solutions at r = R we obtain finally: 
1. In the region where r > R 



e -r I r D 

(r) = Q r D + C . (103) 
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2. 


In the region where 


\ r < 


R 


















<fo(r) = 


Qor 2 D 


+ Co 


e -r /m 

r 


+ 


C 

2ir 


e Y 
R 2 


q \ r 


R- 

e r o 


V 


■0 




with 






Qo = 


Ql 


i + 


Ce Y 
2irR 2 r 2 D 


-A 


R 
r 


R-r 

e r D 



e -r/r D 

Q r 2 D + C • (104) 



lj (105) 

as new non-screened charge arose due the external pressure. The value of Co is determined by 
some additional boundary conditions which we do not consider in the problem. 



The Debye length, Eq. (j97j) . depends on the particles density n and we could define the plasma's 
parameter as 

M = z-7^ < 1 • (106) 

n J TD(n> d 3 x J d 3 pfo(r,p) 

Requiring the validity of this parameter in both cases of presence and absence of external field V2 we 

obtain, that in the case of the existing of external pressure the Debye radius is smaller than in the 

case of the absence of the pressure: 

(r D (Y)) V2 ^ < (r D ) V2=0 . (107) 

It is known, see [llj, that the hydrodynamic description of the process is possible when the charac- 
teristic external parameter of the system, in our case it is R from Eq. ([3]) , is larger than the td '■ 

R > r D . (108) 

Due to the external pressure the Debye radius decreases, and, therefore, we arrive to the inequality 
Eq. (|108p at some large enough value of Y. We see, similarly to the previously obtained results, that 
we can trace the transition from the gas to the fluid state of the system by analytical methods where 
dependence on the external parameters in the conditions of the transition is introduced!!. Therefore, 
perhaps, this approach can be useful in analytical research of the initial stages of evolution of quark- 
gluon plasma at high energy collisions. 

Now, in the linear approximation over fo{r,p) , we have the Boltzmann equation 
dfi(r,p,t) | v dMr,p,t) | / dV 2 (r) | ^ Q \ d h{r,p,t) { g£ll df (r,p) =Q 



dt dr \ dr ' ) dp ' dp 

together with the Maxwell's equations for the electric field: 



rot 
where 



E 1 = 0, divE 1 = Airqn f f 1 (r,p,t)d 3 p, (110) 



E° = -grad(p (r) (111) 

with 4>o(r) given by Eq. (|103|) and Eq. (|104p . Solutions for the fi(r,p,t) distribution function and E 1 
field we are deriving in the Appendix B of the manuscript. 



7 In general, instead rapidity dependence, any other parameters, which characterizes the external pressure and reduce 
the problem to equilibrium one, can be introduced. 
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6 Conclusion 

In our paper we considered a toy model of particles interacting at high energy based on the effect of 
inertial confinement like in plasma physics. There are few main problems which were considered in the 
framework of the model and which can be interesting for the case of quark-gluon plasma description 
as well. These problems are the following: the transition from the gas state to the liquid state of 
the system and a change of the transport properties of the system of interacting particles during this 
transition. 

The trace of the transition to the liquid phase can be characterized by the so called Knudsen 
number in hydrodynamics, see [11] . In our model we demonstrate how much during the evolution of 
the system this number is changing and starting to be largqj , see Eq. (|84p and Eq. (|108p . The largeness 
of this ratio is obtained in two different calculation schemes: by the description of the process in the 
shock waves framework and in the framework of Vlasov's equation. Both approaches are based on the 
assumption that in almost ideal liquid, which is QGP (see [I]), a contribution of dissipative processes 
is small, i.e. dissipation could be manifested on the scales much smaller then the scales considered 
in the paper. In this case the compression of the "internal" gas area by external pressure could be 
considered as an adiabatic process for the convergent to the center of the area of the shock wave. 
This problem of searching for self-similar solutions for the case of non-polytropic gas was considered 
in [8], we adopted the results of this paper for our aims. Also, the assumption on the absence of the 
large-scale dissipative processes allows us to consider the Vlasov's approximation of the Boltzmann 
equation for the self-consistent field for our system as an acceptable as well. The external pressure, 
introduced in the mode in correspondence to the receipts of [3J, reduces the highly non-equilibrium 
process to an equilibrium one in some external field. It is interesting to understand a possibility of 
similar description of the collision process in the framework of QCD, this work we are planning for 
the future. 

As mentioned above, the large value of the ratio Eq. ([84]) and Eq. ()108[) is a sign of the liquid phase 
creation. The dependence of this ratio on the external parameters of the problem is an important 
result, because the main tool in quark-gluon plasma evolution is the hydro dynamical description of 
the plasma state, see [IJ. Therefore, our approach in application to QGP can serve as a description 
of the initial states of the hydrodynamical evolution. Namely, it will be interesting task to use the 
proposed approach for the definition of the initial states of hydrodynamical expansion in the framework 
of Landau model, see [I] and for calculation of multiplicities of the produced partons in high energy 
interactions, see [32] and references therein. 

The QCD QGP is an almost ideal fluid, where interactions with large coupling constant between 
the participants of collision are revealed, see [1] for review. This fluid has very special transport 



8 Our ratio is opposite to the usual Knudsen number, therefore in our approach the liquid state is characterized by 
the large value of the ratio. 
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properties, which so far can be partially explained only in framework of AdS/CFT duality model; 
only there the calculations in the strong coupling constant domain are possible. Nevertheless, in our 
approach we tried to consider the subject of the transport properties of the system under the external 
pressure using well known Enskog calculation scheme, see Eq. (|2Up and [7j. It is important to underline 
also, that in Enskog approach we calculate the corrections to the transport coefficients caused by the 
external influence on the gas whereas the initial values of the parameters are given. Our results 
demonstrate that in this case the transport properties of the system of interacting particles near the 
transition into the liquid phase, see Eq. ([Mjl -Eq. (|36|) . are similar in some sense to the properties of the 
quark-gluon plasma obtained in [S]. Namely, we obtain similar rise of bulk viscosity when temperature 
is rising. The explanation of this effect in our model is simple. In Enskog approximation the bulk 
viscosity is proportional to the density of the state. The density is growing toward some constant 
value whereas the system of interest undergo a gas-liquid transition. When the state of ideal liquid 
is created finally then the density is constant due the incompressibility of the state and, therefore, 
the bulk viscosity as well. Another interesting result is an equation Eq. (|39|) ; we see that this ratio is 
small in our framework similarly to what is observed in QGP , see for review and explanations pQ. 

Finally we conclude, that the main purpose of the exploration of our toy model was an under- 
standing of the main principles and methods of the approach for it's generalization at the case of 
quark-gluon plasma. There are the following discussed problems which we found also interesting to 
investigate at the case of the quark-gluon plasma: 

a The dynamical change of the Debye length of the plasma under the external pressure. This 
change indicates the gas-liquid transition by the change of the Knudsen number of the system. 
In our paper the Debye length change dynamics was depend on the model introduced potential 
for the external field pressure. Nevertheless, if we have a plasma state under some external 
pressure the similar effect must be observed for any kind of the plasmas. 

b Calculations of the one-particle distribution function and self-consistent fields in the Vlasov's 
approximation for the dense compressed drop of the partons in the field of external particles. 
This question is especially interesting because this function can clarify the transport properties 
of the system which are very unusual at the case of QGP. 

The obtained results show that these tasks can be achieved in the approach correctly reformulated for 
the case of QCD QGP that we plan as our future work. 
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Appendix A: gas of interacting particles 

In this appendix we shortly remind the main facts concerning a virial expansion, see the detailed 
derivation in [T3] . We consider a gas of N interacting particles each with mass m on a plane as a gas 
of hard disk with small thickness and radius tq. The energy of the gas in the classical limit is given 
by the well known expression: 

N 2 

E(p,q) = J2^ + U(h,...,b N ) (A.l) 

where as usual the first term is the kinetic energy of N particles, U is a potential energy of their mutual 
interactions and b\, ..., 6jv their coordinates. The grand partition function for this Hamiltonian is 

oo -, „ N t2 - j2T N 2 

Q(»,T,S) = e -/»" = £ ^/n ^exp(-^|) exp(-/3 t/(6 1; ...,&„)) . (A.2) 

7V=0 ' i=l i=l 

where is the grand potential of the problem. As usual, integrating over momenta, we obtain 

Q(,,T,S) = e -/>° = £ j^)" M^> (A.3) 

and 

Z N (S, f3) = J ... J d%, ..., d 2 b N exp (-/3 [/(fox, ..., 6jv)) . (A.4) 

Here A is the De Broglie length of the quark corresponding to the average energy f3: 

A=(-^-f. (A.5) 

\2irmJ v ' 

The potential of the problem is given by 

n = -^ln(l + ^S + ^ J J d 2 b 1 d%ex V (-/3U(b 1 ,b 2 )) +...J . (A.6) 

Here, we used 

i2u o d2 



dfb = S = vr R z , (A.7) 

where R 2 is the characteristic radius of the problem. In the following we will define and consider only 
pairwise interaction between the particles, namely we have 

tf (61,62) = U(\h -6 2 |) = C/(6i 2 ) = C/(6) = *7i2 (A.8) 

and therefore, in the relative coordinates of the center mass we reduce the multiplicity of the integrated 
functions and obtain an additional S factor in the integrals: 

n = -P5 =-- In f 1 + S—+ <?^ J J d 2 b 12 exp(-pU 12 ) + ... J . (A.9) 
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Introducing variable £ 



C = -jp (A.10) 



we obtain the expression for the potential in the form of the series in £ 

n = -ps=-f E^r- (A.11) 

^ n=l 
We will take into account only two first terms of this series with the following J\ and J<i'- 

Ji = 1, J2= I I d% 2 (exp(-/3 17 12 ) - 1). (A.12) 



The number of particles in this gas we obtain as usual 



\dnJ T>s 

and because d( / d fi = f3 (, we finally have: 

OO x 

Excluding from Eq. (jA.llj) and Eq. (JA.14J) the variable (, we obtain in the second order approximation 

the equation of state for our gas: 

NT N 2 T 



P = — ~YW J ^ (A ' 15) 



where T = 1//3. 
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Appendix B: deviation from equilibrium state and non static electri- 
cal field of the gas in Vlasov's approximation 

Proceeding with the equation Eq. ()90p we perform a following substitution: 

fi(r,p,t) = fi(r,p,t) f (r,p) (B.l) 

with fo(r,p) from Eq. ()93j) . Rewriting Eq. f)109|) we obtain: 

d fi(r,p,t) p dfi(r,p,t) „dfi(r, p, t) a E 1 . , 

a l + - ■ /U ^' ; +F ^'^ ; - * ^ P = 0, (B.2) 

ot m or op mkB ± 

where F is given by Eq. (|87|) . We consider a linear approximation over the equilibrium distribution, 
therefore, we could write the external force in the equation in following form 

F(r) -»■ F = F(f) = - ( dV ^ ) +qE°(r). (B.3) 

with r 

_ = _ / rf (r,p)d 3 xd 3 p 

I fo(r,p)d 3 xd 3 p 
Thereby we have the following equation for the fi(r,p,t) distribution function: 

dfi(r,p,t) p 9/i(r,p,t) -d/i(r,p,i) g^ 1 

«7 ' « h F « j — ^r P = ° ( B - 5 ) 

at m ar op mkB± 

Performing Fourier transform of fi(r,p,t) and E l (r,t) 

fi(r,P,t) = J^-l-^e-^ t+irk Mk,P^) (B.6) 

we obtain finally a nonlinear differential equation of the first order over momenta: 

d<t>o 



dp 
Here 



+ apcpo — b(j)Q — cp = . (B.8) 



ik , iuj a E 1 , n _ 

= — F ' 6 = ~P ' C = I TF • (B - 9) 

7TT, i< i* mkB 1 b 

The solution of this equation, with additional condition that at p — > (T —} ) the solution is static, 
i.e. (fto(k,p = 0,u>) = , is the following function: 

^o(Jfe,p,w) = ce - a P 2 /2 + b P [ P p e ap ' 2/2 - bp ' dp . (B.10) 

Jo 

Thereby we obtain for the correction to the equilibrium state: 

S-^ 1 _-ap 2 /2 + bp f „' ap'2/2-fep' j ' 



0o(&,^) = | -^ e-" /^ / p e ap >*->P dp . (B.ll) 

rnksT F Jo 
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We see, that this correction is suppressed in comparison to the distribution function Eq. (|93p by the 
factor 2-gr-. Substituting this expression back in Eq. (|B.6p . we have: 

(B.12) 
Performing Fourier transform again, we obtain for our function in (r, t) representation: 

fi(r,p,t)= q =, ! dp'pE 1 (r + P —JL,t + P -^). (B.13) 

m Kb J- r Jo 2m b b 

This correction to the static distribution function determines also the non-static electric field which 
we consider further. 

Using Eq. (|89p . we find the equation for the correction E 1 to the electric field: 

divE\r,t) = 4 ^ 2 " / d 3 p f dp'pE\r + ^—^,t + tJZJL) f^ p ) (B.14) 

rnks 1 r J Jo 2 m b b 

with with fo(r,p) from Eq. (|93p . see also Eq. (|B.lj) . This equation is highly non-linear and non-local 
and, perhaps, a precise solution of the equation can be found only by numerical methods. In order to 
investigate approximate solutions of the equation, we perform Fourier transform of the functions in 
Eq. tiB~T4l) : 

,*l„ x Wn f j3 fP'>fd 3 k'~ ln , --feW- fc ')(4#) ~ n , , 
kE\k,u) = —^g I d 3 p J dpp J 7 —-E 1 {k-k,u})e V ) \ Jf (k,p), 



mk B TF J " Jo J (2tt) j 

(B.15) 

where 

d 3 k 



fo(r,p) ~- J j^e" k fo{k,p). (B.16) 

In the right hand side of the equation the oscillating integral over k is not vanishing when 



2mF F 

p" 2 — p 2 k B T 



k oc 2 ex ^^-^ « 1 (B.17) 



in our approximation of the weak external field. Therefore, in the region where 

k > k' (B.18) 

in the first approximation over k we have: 



kE 1 (k,u})faE 1 (k,u) mk H T - j d A p J dp p J -—e \ J V //o(*,p)- 

-^hr^^fEj dp Jo dpp Jj^ ke v y v y/o(jfe ' p) (R19) 



22 



or, simplifying the notations we obtain: 



dEl Q k k U) ~ E\k,u) (C_(w,fc) - C 2 (u,k)k) = 0. 



(B.20) 



where 



»__; \ +ik l -— i 2 



JfpJSdp'p'e VW V— ' M*^r,P) 



„2 J 2 



ciK^) = — u 7-— — , ;;,: 2 rc — ( b - 21 ) 

f<PpJSdp'p'f-$frk'e 



n> ,/„'„' , •^z.v^^)^*-*'^^)/,,^) 



and 



C2(«, *) = m^TF _____ . (B22) 

4 7 r^ n/ ^ /(f ^ V/ |fe; fc ' e --l £ ^J +l(fe - fc) l" )f (k',p) 

An integration of Eq. (|B.19j) gives 

E\k,u) = C (u)efi Cllu ' t)dt -ti Oa{u ' t)tdt . (B.23) 

We see, that the fluctuation of electric field in this region of k is suppressed by the large factor C2 in 
the power of exponent. 



In the opposite limit, when 



, 2mF F 

k ~ k oc -75 ^ oc — — < 1 (B.24) 



a l \ 1 ;. i p 2 ~p 2 



we introduce new variable e: 

e = k — k . 

Equation Eq. (|B.19P will acquire the following form: 

Fourier transform of function Eq. (|9_j) is 

f q ■ fc _Z2____) / „ _ V 9 (k) Q(bn(k)\ 

/o(M = /„(*)/o(p) = /o(p)y Ae- 4rfc e *_* *_* « / (p) f (2^) 3 5 3 (A ; ) - -j±± - ^^j 

(B.26) 

with 

1 



fM ~ ^^¥W e • {B27) 



therefore, in the first approximation of expansion over e we obtain: 

E\k,u) = E\k,u) t'C/L f d3 Pfo(p) [ P dp'pe 
rnks 1 r k J Jo 



W™ r , 3 ,, , /* , , - fa l E f+«M^ 



2 „2 
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-i, x 47r« z nK 2 fc f ,, „ , , f p , > > ~ lw \ f^ f 



m(k B T) 2 FkJ " JVK " J ** J 2vr 2 



4vr<f n0 o (/c) /" )S , , N T , ' ' ~ iw l V* ) /"fcfr dee 2 .„ , 



Finally we obtain for our electric field: 



~i„ ^ ~i/ s 2q 2 nF 2 ( V 2 (k) qMk)\ f ,, „ /x P , / / "^(vO 



where 



(B.29) 



£ " * - ^TFk I dpfM I """ e ' ' V ' (B30) 



4vr q 2 n f 3 £ f _ s [P _,_, _, - iw [ V^ j +ik \ ■•„, , 



£ L l_i_, (. I P 2 -P 2 



is the dielectric constant of the problem. 
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